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I.  INTRODUCTION 


With  the  development  of  LDV,  the  flow  of  dilute  gas-solid 
suspension  has  become  of  considerable  importance.  A  simple 
situation  of  such  a  flew  is  the  motion  of  a  single  {.article 
in  the  turbulent  stream.  The  most  important  question  which 
arises  in  such  a  situation  is  whether  the  solid  particle 
follows  the  streamlines  or  not.  Since  LDV  measures  the  shift 
from  a  single  scatterer  it  is  of  importance  to  evaluate  as  to 
how  far  the  results  of  velocity  as  well  as  that  of  frequency 
spectra  of  tracer  particles  provide  the  requisite  information 

regarding  the  turbulent  structure  of  the  gas. 

1  2 

Tchen  (see  Ref.  3),  Hinze  ,  Corrsin  and  Lumely  ,  and 
3 

Soo  have  examined  this  problem.  Starting  with  Tchen 's  equation, 

4 

Hjelmfelt  and  Mockros  have  studied  the  importance  of  added 
mass,  history  of  acceleration  and  pressure  gradient,  due  to 
fluid  acceleration  in  the  motion  of  a  particle.  The  analysis, 
based  on  a  formulation  due  to  Tchen,  led  to  the  conclusion 
that  the  particle  diffusivity  is  the  same  as  the  Lagrangian 
eddy  diffusivity  of  turbulence.  Later  experiments  by  Soo 
and  studies  based  on  'probability  of  encounter'  show  that  these 
two  dif fusivities  are  different  and  the  particles  do  not 
follow  the  fluid  particles.  Even  for  a  very  small  particle 


1 


which  is  expected  to  follow  the  flow,  the  increase  in  flow 
Reynolds  number  decreases  the  particle  diffusivity. 

f 

In  the  present  report,  the  theory  of  Tchen  is  extended  to 
incorporate  the  difference  in  these  two  dif fusivities,  as  well 
as  the  effect  of  temperature,  density,  and  particle  size  distri¬ 
bution  on  LDV\measurements , 

*>» 

II.  BASIC  EQUATION 


Based  on  the  assumption,  that  the  presence  of  a  particle 
does  not  modify  the  flow  and  that  the  relative  motion  gives 
rise  to  Stokes  drag  on  the  particle,  Tchen  formulated  the 
following  equation  for  the  motion  of  the  particle: 


du 


dt 


=  ~^£- 


2  p  +  p 

PP  P 


du . 
r — i. 
“dt 


2 

3 


V 


ax .  ax 

•J  J 


+  A  ^  (ui-v+  p<vv  -&] 


+  -—2 — ~  f  P  dT(ui~up-i  ) 

(2p  +P)a  v  \  J  - - cJt 


(1) 


tfhere  u  is  the  vector  velocity  of  the  particle,  u.  is  the 
i  i 

velocity  of  the  gas,  v  the  coefficient  of  kinematic  viscosity 
u  the  coefficient  of  viscosity,  p  the  density  of  the  gas,  p 

P 

the  density  of  the  solid  material  and  a  is  the  radius  of  the 


2 


particle.  Tchen  assumes  that  the  particle  radius  is  small  and 


u 


»  i 


(2) 


The  most  restrictive  assumption  mady  by  Tchen  is  that  during 
the  motion  of  the  particle,  the  same  fluid  element  remains 
in  its  neighborhood.  Thus  he  postulates 


d _ rd_ 

dt  dt 
P 


The  consequence  of  this  assumption  is  that  the  streamlines 
and  the  trajectory  of  the  particles  coincide  and  as  such, 
leads  to  the  same  diffusivities  for  the  particles  and  gas. 
In  order  to  account  for  the  decrease  in  diffusivity  of  the 
particles,  we  introduce  a  relaxation  process  such  that 


dUpi  du 
dtp  -  dt 


+  vu 

P 


(3) 


where  y  is  determined  such  that  the  appropriate  dif.  usivity 
is  obtained  in  the  limit  of  large  t.  Utilizing  (2)  and  (3) 


in  Eq.  (1  ) ,  we  get 
du 

d^  +  S  dt 

where 


+  vu  =  B  ^  +  aB(u-u  l  +  Bf^  f 


— (u-u  ) 

P- 


(4) 


-ao  Jt-T 


a  =  3v/as  and  0  =  3  p/ ( 2  p  +  p ) 


Following  Chao  (see  Ref.  3),  „e  define  the  Fourier  transform 
of  any  quantity  as: 

poo 

“  I  9(t)exp(-iuut)dt 

v -oo 

Where  for  non-decaying  turbulence  the  above  integral  is  taken 
to  be  in  the  sense  of  generalized  harmonic  analysis.  Taking 
the  Fourier  transform  of  Eq.  (4)  leads  to 

un  (i+/Ti)+iu  rr^l 

ZR  -  -  -  v  2  a'  La  +  V  2  a  1 _ 

u  " '  ' 

(l+  +  /I  Hi)  +  ,-/“L  4.  /~\ 

v  a3  V  2  a'  +  M  ag  +  V  2  a ) 


where 


.  £ 

u3 


■(wt?)  *(5vT!) 

■  f1  ♦  VVH)’ ♦  fe  ♦,/!« )' 
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It  may  be  noticed  that  as  ou/a-*0,  _  - 1-  - 

'  n3  1  +  Y/aP  and  UJ/ot-,°° * 
Q , 

■  1  — »  Q  2 
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III.  AUTO  CORRELATION 


The  Lagrangian  velocity  auto  correlation  for  the  gas 

i 

is  given  as 


R  (  T )  =  <U(t)u(t+T)> 

<U3> 


where 


T 

<u(t)u(t+T)>  =  lim  ±-  [  U(t)u(t+T)dt 
T^oo  J  m 


-T 


The  spectrum  density  F(uu)  is  defined  as: 
<u*>F(uu)  =  lim 


V T 


where 


T-<x> 

— i  uut 


2TTT 


=  J  u(t)e"llljtdt 


(6) 


-T 


and  star  denotes  the  complex  conjugate.  The  Wiener -Khintchine 
theorem  gives  a  relationship  between  R(t)  and  F(ou)  as 


R(t)  =  J°° F  ( uu)  cos uut  duu 
o 


woo 


F  ( a))  =  —  J  R{t)cosajt  dt 
o 


(7) 


A  similar  set  of  relationships  can  be  introduced  for  the 
particles.  These  are  denoted  by  a  subscript  'p'.  Thus 


<ua> 


f 


lim  ( 


u  u 
P.T  p, T 


o  T-oo 


2  TTT 


using  Eqs.  (5)  and  (6),  we  get 


2  S  O 

<u  >  =  <u  >  I  F  ( U))  duj 


The  spectrum  density  function  for  the  particles  is  then 
given  by 


V”*  =  o'  p<»'> 


(8) 


Equation  (8)  clearly  shows  the  relationship  between  the 
spectral  density  measured  by  LDV  and  the  spectrum  of  turbulence 


of  the  gas.  It  may  be  pointed  out  that  the  above  relation  is 
valid  for  particles  of  single  diameter.  The  presence  of 
various  size  particles  in  the  flow  further  complicates  the 
simple  relationship.  If  f^a)  is  the  distribution  function 
for  particle  size,  then  Eq.  (8)  is  replaced  by 

V'">  '  F(») £,,(»)<!» 

where  ax  and  a#  are  the  smallest  and  largest  particle  radii 

present  in  the  flow.  It  is  the  'id'  dependence  of  and  the 

‘‘a 

radius  ’a'  which  limits  the  usefulness  of  LDV. 


IV.  DIFFUSIVITY 


For  homogeneous  turbulence,  it  has  been  shown  by 
Kampe'de  Feriet 

D(t)  =  <u3>  j  R ( t) dt 


Similarly, 

Dp(t)  =  <u>  J  VT)dT 
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utilizing  Eq .  (7)  and  interchanging  the  order  of 
integration,  we  get 


Vt,  =<u;>f  s4a*pp,„d. 


Thus 


sinuut  _  .  .  . 

l£  =  _Ze.  .  'o  ~T~  Fr>»dui 

D  <U3>  jx 


<u 


1“  Siiiat  F(lIl)dlll 


(9) 


«j  uu 

o 


Substituting  from  Eq.  (8)  in  Eq.  (9)  and  taking  the  limit 


t-<co  ,  we  get 


.  .  D 

lim  _jd  _  1 

t-oo  D  1+y/ag 

For  finite  y ,  the  last  equation  exhibits  the  desired  reduction 

m  the  diffusivity  shown  by  experiments.  As  shown  in  Eq .  (5), 

it  may  be  seen  that  for  a}=0, 
u3 

p  _  1 

u  i+y7ob 

i.e. , the  two  speeds  are  not  the  same.  Since  at  large  Reynolds 
number  the  particle  diffusivity  decreases  (even  for  small 
particles) ,  the  particles  in  such  a  flow  may  not  follow  the 


streamlines  at  all. 


V.  RESULTS 


In  Figures  1-4  the  expression  (5)  for  up2/u2  is  plotted 
against  u)  for  various  s’ze  particles  and  diffusivity  ratios 
shown  on  the  graphs.  In  the  limit  of  a-o,  y-'oo ,  B^o  with 
ag-const,  Eq.  (5)  reduces  to 


l+(- 


9  4 


3 


which  is  the  expression  for  the  motion  of  a  particle  acted 


upon  by  Stokes  drag  only.  The  values  of  Pp*P  and  v  are 
taken  to  be  2.25,  .00118  and  . 157, respectively .  If  we  introduce 
the  Cunningham  correction  into  the  coefficient  of  viscosity, 
the  above  expression  may  be  written  as 


u' 


JE  = 


u' 


l+t 


r2TTfp  Ds 
f  -  P  P 


184(1+  £i) 

Dp 


I 


(10) 


where  D^-2a,  UJ=2rrf  and  1  is  the  molecular  mean  free  path  and 
k  is  the  Cunningham  constant  which  is  1.8  for  air.  Expression 
(10)  has  been  used  by  Becker,  Hottel  and  Williams5  and  many 
other  authors.  In  order  to  assess  the  effect  of  variations 
in  density  and  temperature  of  the  medium,  which  enters  expres¬ 
sion  (10)  through  4  and  t,  isentropic  relationships  are  used 
and  these  variations  are  parameterized  by  introducing  the  Mach 
number  and  stagnation  temperature.  Figure  5  shows  the  plot 


8 


of  up^u  as  a  Action  of  £  for  a  Mach  number  of  0.26  and 
various  size  particles.  The  stagnation  temperature  for  this 

i 

case  is  taken  to  be  the  atmospheric  temperature.  Figure  6  is 
a  similar  plot  for  Mach  number  10  and  the  stagnation  temperature 
ratio  of  about  4.7.  In  Figure  7,  the  effect  of  variations  in 
stagnation  temperature  are  examined  for  a  flow  of  Mact  number  2. 

In  order  to  investigate  the  effect  of  the  presence  of 
various  size  particles,  a  distribution  function  of  the  following 
form  is  used: 


£«<a)  =  A*exp-S,s 


where  Ax_  is  chosen  such  that 
N 


da  =  1 


and  AA  is  the  measure  of  the  flatness  of  the  distribution 


u 

function.  The  integral  of  f  (a)  over  'a'  for  a  =1.0U, 

U  IN  O  M 

Dr=  1.0,  ai=0.5a  and  as=3.5^i,  and  for  various  values  of  AA 

(shown  on  the  figures)  is  shown  in  Figure  8.  In  Figure  9,  the 
r  u 

argument  of  J  fN(a)  as  a  function  of  ui  is  shown.  The  distri¬ 
bution  function  is  plotted  in  Figures  10  and  11.  From  Figures  8  and 
9,  the  fall  in  the  response  of  the  particles  is  quite  clear.  It  may  be 
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mentioned  that  DR,  the  ratio  of  the  diffusivities  is  also  a 
function  of  the  particle  size.  Although  a  relationship 
giving  the  particle  size  dependence  and  Reynolds  number  is 
yet  to  be  established,  but  its  influence  on  the  particle 
motion  when  various  size  particles  are  present  can  be  clearly 
deduced  from  Figures  1-11. 

VI.  CONCLUSION 

Motion  of  a  particle  in  a  turbulent  flow  has  been 
examined.  It  is  found  that  the  frequency  spectral  density  of 
th*_  particle  motion  depends  upon  many  factors;  notably  the 
particle  size,  ratio  of  particle  to  fluid  diffusi  'ity  and 
the  Mach  number  of  the  flow.  Increase  in  particle  size  and/or 
Mach  number  and  decrease  in  particle  diffusivities  leads  to 
a  smaller  value  of  frequency  of  turbulence  to  which  particle 
will  follow  the  flow  field.  This  fact  restricts  the  applica¬ 
bility  of  LDV  to  low  frequency  and  low  .speed  turbulent  flows, 
assuming  that  a  single  particle  is  present  in  the  scattering 
volume.  The  uncertainty  in  the  size  distribution  of  the 
naturally  present  particulates  in  a  flow,  coupled  with  the 
uncertainty  in  the  distribution  of  the  artificially  introduced 
scatterers,  their  coagulation  properties,  their  diffusivities. 


10 


etc.,  limit  the  field  of  applicability  of  the  LDV,  in  particular 

to  turbulent  measurements.  This  suggests  that  for  a  more 

« 

meaningful  interpretation  of  the  LDV  measurements,  a  monitoring 

of  the  sizes  of  the  particulates  should  be  carried  out  and 

* 

incorporated  in  the  data  reduction  process. 

A  preliminary  assessment  of  the  presence  of  high  tempera¬ 
tures  and  changes  in  density  and  temperatures  has  been  presented 
in  Figures  5,  6,  and  7.  Radiometric  forces  due  to  the  presence 
of  temperature  gradients  and  radiation  transfer  should  also 
be  included  in  Eq .  (1).  The  later  forces  couple  the  motion 

of  the  particle  to  the  temperature  fluctuations  in  the  gas. 
Although  these  are  higher  order  effects,  they  may  further  limit 
the  applicability  of  the  LDV. 
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FIG. I  RATIO  OF  PARTICLE  ANO  FLUID  VELOCITY  AS  A  FUNCTION  OF  FREQUENCY 
FOR  VARIOUS  VALUES  OF  DIFFUSIVITY  RATIO  FOR  a-0.000025  cm. 
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FIG.  8  EFFECT  OF  THE  PRESENCE  OF  VARIOUS  SIZE  PARTICLES  ON  THE 
RATIO  OF  PARTICLE  AND  FLUID  VELOCITY  AS  A  FUNCTION  OF 
FREQUENCY 


